Normalitity preserving perturbations and augmentations and 

their effect on the eigenvalues 



Ricardo Reis da Silva and Jan H. Brandts 

O 

psj January 13, 2013 

Abstract 

We revisit the normality preserving augmentation of normal matrices studied by 
I— I Ikramov and Eisner in 1998 and complement their results by showing how the eigen- 

values of the original matrix are perturbed by the augmentation. Moreover, we construct 
all augmentations that result in normal matrices with eigenvalues on a quadratic curve in 
the complex plane, using the stratification of normal matrices presented by Huhtanen in 
2001. To make this construction feasible, but also for its own sake, we study normality 
preserving normal perturbations of normal matrices. For 2x2 and for rank-1 matrices, 
^ the analysis is complete. For higher rank, all essentially Hermitian normality perturba- 

' tions are described. In all cases, the effect of the perturbation on the eigenvalues of the 

^ original matrix is given. The paper is concluded with a number of explicit examples that 

^ illustrate the results and constructions. 

I Keywords: structured perturbation theory, normal matrix, essentially Hermitian matrix, 

Toeplitz decomposition, normal augmentation problem. 

cn 

O 1 Introduction 

. . A well-known sufficient condition for the sum of normal matrices A and E to be normal is, 

that A and E commute. This condition, however, is not necessary. If A and E are both 
^ Hermitian, the sum is also trivially normal. In both the above situations, a well-developed 

perturbation theory for the eigenvalues oi A + E in terms of A and its perturbation E is 
available. See for instance p2 j and the references therein. For other normality preserving 
normal perturbations, one could of course apply one of the many equivalent characterizations 
[21 |6j of the normality of A + E, but other than that there seems to be much less literature. 
Although Wielandt [13] already studied the location of eigenvalues of sums of normal matrices 
in 1955, he did not require the sum to be normal. 

In this paper, we study the perturbation of normal matrices by essentially Hermitian 
matrices [HEKH]. These are matrices E that can be written as E = f3H + al, where q,/3 G C, 
H is Hermitian, and I is the identity matrix. Since all 2 x 2 normal matrices and all rank-1 
normal matrices are essentially Hermitian, this will completely characterize the 2x2 and the 
rank-1 case, but will also give insight in the k x k and rank-(A; — 1) case for A; > 3. Moreover, 
we will relate the eigenvalues of ^ -|- to those of A. 

Essentially Hermitian matrices have been frequently discussed, both in the numerical and 
the core linear algebra setting. See for instance the introductory section of [5], where the dis- 
cretization of the Helmholz equation Au + k'^u = f where k = ki-\-k2i € C is the wave number 
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and u the unknown eigenmode of the differential operator, results in essentially Hermitian 
matrices. Faber and Manteuffel show in |4j that for linear systems with essentially Hermitian 
system matrix, there exists a variant of the Conjugate Gradient method that still relies on 
a three term recursion. In the context of the Arnoldi method applied to normal matrices, 
Huckle [7J proves that irreducible normal tridiagonal matrices are essentially Hermitian. In 
fact, the upper Hessenberg matrices generated by the Arnoldi method is tridiagonal only for 
essentially Hermitian matrices. The behavior of Arnoldi's method is in that case similar as 
the Lanczos method. Generally, however, the Arnoldi method for normal matrices does not 
respect the structure. Huhtanen in |H] raises this issue. He shows that for every normal matrix 
A and almost all unimodular z £ C the skew-Hermitian part S of zA is a polynomial p of 
degree at most n — 1 in the Hermitian part H of zA. The essentially Hermitian matrices are 
exactly the ones for which the degree of p is at most one, resulting in collinear eigenvalues. 
Since the p can be retrieved in a modest number of arithmetic operations as a by-product of 
the Lanczos algorithm applied to H, this led to the development [HI [9] of efficient structure 
preserving algorithms for eigenvalue problems and linear systems involving normal matrices. 

Essentially Hermitian matrices also appear in the context of core linear algebra problems. 
For instance, in [IQ], Ikramov investigates which matrices R can be the upper triangular part 
of a normal matrix A. Although the full problem remains open, a solution exists if all diagonal 
entries of R lie on the same line in the complex plane, say 

^:R^C: p^y + Op, y G C and 9 e C, \9\ = I. (1) 

Writing R = D + U where D is the diagonal of i?, we see that 9[D + U — yl) is upper triangular 
with real diagonal, and consequently, H = 9{D + U — yl) + 9U* is Hermitian. Therefore, 

A = 9H + yl = R + 9'^U* (2) 

is essentially Hermitian, hence normal, with upper triangular part equal to R. Together with 
Eisner [llj, Ikramov then considered the normality preserving augmentation problem: when 
does the addition rows and columns to a given normal matrix A yield a normal matrix A^? 
Essentially Hermitian matrices appear in this setting as special cases, though the problem in 
its full generality remains unsolved. We will take the augmentation problem as a starting point 
for this paper. We add to their results an analysis of the eigenvalues of the augmented matrix 
yl+ in terms of A. Moreover, we investigate how Huhtanen's stratification of the normal 
matrices may be of help in providing additional solutions. Our analysis involves normality 
preserving normal perturbations as mentioned in the beginning of this introduction, a topic 
that we will study in more detail as well. Also in this context we present results on the 
eigenvalues of the perturned matrix. 

1.1 Outline 

In Section [2] we define counterparts of the Hermitian and skew-Hermitian parts of a matrix, 
based on unimodular complex numbers 9 and 9i, that for 9 = 1 reduce to the standard 
Hermitian and skew-Hermitian part. This will somewhat facilitate the discussion on essen- 
tially Hermitian matrices of the form 9H + al. We also outline the description of normal 
matrices as given in [8j and pay attention to commuting matrices and their eigenvalues and 
eigenvectors. In Section [3j we revisit the normality preserving augmentation problem by 
Ikramov and Eisner [3] and add to their analysis results on the eigenvalues of the resulting 
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augmented matrix. In Section [4] we use similar ingredients to study the normahty preserving 
normal perturbation problem. This will be of use in Section [5| where we explicitly construct 
augmentations of normal matrices that have all eigenvalues on a quadratic curve in C. 

1.2 Notation 

For quick reference, we list here some notations that are frequently used in this paper. 

• M"^™'(C) is the set of all n x m matrices with complex entries, 

• /, or /„ G M"^"'(C) is the identity matrix with columns ei, . . . , e^, 

• cr{A) the set of eigenvalues of A, its spectrum, 

• T is the circle group of unimodular numbers, 

• Tu is the subset of T in the closed upper half plane, minus —1, 



• @{A) and @ (A) are 0-(skew)-Hermitian parts of A (see Section 2.1), 

• [A, B] = AB — BA is the matrix commutator. 

• V^{M.) is the space of real- valued polynomials on M of degree less than or equal to k 

The zero matrix we simply denote by regardless of its dimensions, and within matrices we 
will often use empty space to indicate a zero block. 

2 Preliminaries 

In this section we define a family of non-standard decompositions of complex numbers and 
matrices and prove some of their properties. We also recall a result on the eigenvalues and 
eigenvectors of commuting normal matrices. 

2.1 The ^-Toeplitz decomposition and some relevant properties 

A complex number z G C is often split up as z = + where 25}?(z) = z + ^ and 

29(z) = z — z. This interprets the complex plane as a two-dimensional real vector space with 
as basis the numbers 1 and i. In this paper, it will be convenient to decompose z differently. 
For this, we introduce a family of decompositions parametrized in ^ G T C C, the circle group 
of unimodular numbers. For a given G T we let 

z = e(z) + e^(z), where Q{z) = ^(ez)9 and Q^{z) = Q(ez)9. (3) 

Moreover, apart from the standard Toeplitz or Cartesian decomposition of a square matrix 
A into its Hermitian and skew-Hermitian parts, 

A = 'H{A)+S{A), where n{A) = '^{A + A*) and S{A) = ^{A - A*), (4) 

in accordance with ([s]), we consider the family of matrix decompositions 

A = Q{A) + Q^{A), where e{A) = n(9A)e and O^iA) = S(9A)9. (5) 
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We will call this decomposition of A its 9-Toeplitz decomposition, Q(A) the 6-Hermitian part 
and Q-^{A) the 6 -skew-Hermitian part of A and subsequently, A is 6-Hermitian if ^ = 0(^4) 
and 6 -skew-Hermitian if ^ = Q-^{A). Note that 9-skew Hermitian matrices are i^-Hermitian. 
We now generalize a well-known result to the 0-Toeplitz decomposition of normal matrices. 

Lemma 2.1 Let 9 & T be arbitrarily given. For any normal matrix A we have that 

Av = Xv ^ e{A)v = e{X)v and 0^ {A)v = {X)v , (6) 

where Q and 0"*" and their relation to 9 are defined in 

Proof. Let Au = Xu for some u ^ 0. Since A is normal, there exists a unitary matrix 
U with u as first column and U* AU = A diagonal. But then U*A*U = A*, showing that 
A*u = Xu. Since the argument can be repeated with A* instead of A, this yields that A 
and A* have the same eigenvectors and that the corresponding eigenvalues are each other's 
complex conjugates. Therefore, 

2e{A)v = (OA + 9A*)9v = (9X + 9X)9v = 2Q{X)v (7) 

and, similarly, also @-^{A)v = Q^{X)v. The reverse implication is trivial. □ 

Corollary 2.2 Let Xi be an eigenvalue of a normal matrix A. For given 9 gT, consider the 
line i C C through Ai defined by 

i -.R^C : Xi + p9. (8) 

Assume that Xi,...,Xp are all eigenvalues of A that lie on i. Then the eigenspace U of 
the eigenvalue ©"""(Ai) of Q^{A) equals the invariant subspace of A spanned by ui, . . . ,Up. 
Restricted to U, the matrix A — XiLp is 9-Hermitian. 

Proof. Let Xa, Xb be eigenvalues of A, then we have that 

e^iXa) = G^(Afc) ^ 9(^Aa) = Qi9Xb) ^ e{Xa -Xb) = p£R ^ Xa-Xb = 9p. (9) 

Thus, @-^{A)uj = @-^{Xi)uj for all j £ {I, . . . and conversely, if Q-^{A)u = Q^{Xi)u then 
n is a linear combination of ui, . . . , Up. Writing Up for the matrix with columns ui, . . . , Up we 
moreover find that 

AUp = UpAp with Ap = 9R + XiLp (10) 

where Ap is the px p diagonal matrix whose eigenvalues are Ai, . . . , Ap and R is real diagonal. 
Thus, Q-^{U*AUp — XiLp) = 0, proving the last statement. □ 

2.2 Normal matrices with all eigenvalues on polynomial curves 

If A is 0-Hermitian then A is normal because 9A is Hermitian, and by the spectral theorem 
for Hermitian matrices, all eigenvalues of A lie on the line ^ : M — )• C : p i— )• p9. In the 
literature, for instance in [U [5] , ^ is called essentially Hermitian if there exists an a € C 
such that A — aL is 0-Hermitian for some 9 £ T. Clearly, the spectrum of an essentially 
Hermitian matrix lies on an affine line shifted over a £ C. Conversely, if a normal matrix 
has all its eigenvalues on a line ^ C C, it is essentially Hermitian. This includes all normal 
2x2 matrices and all normal rank one perturbations of al for a S C. Larger and higher 
rank normal matrices have their eigenvalues on a polynomial curve C C C of higher degree. 
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2.2.1 Polynomial curves of degree k>2 

Each normal A G M"^"'(C) has its eigenvalues on polynomial curve C C C of degree k <n — \. 
This can be explained as follows [8]. First, fix € T such that for each pair Ap, \q of eigenvalues 
oiA 

Ap / A, ^ e(Ap) / e(A,). (11) 

Note that there exist at most n{n + 1) values of for which this cannot be realized, corre- 
sponding to the at most \n{n + 1) straight lines through each pair of distinct eigenvalues of 



A. Once (11) is satisfied, the points 

(^e(Ai),i^e^(Ai)), . . . , {m{\n),iee^{Xn)) e m x m (12) 

form a feasible set of points in M x M through which a Lagrange interpolation polynomial 
TT G P"'~^(]R) can be constructed that satisfies 

e^(Aj) = i-07r(^G(Aj)) for allj G {l,...,n}. (13) 

Notice that 6Tr{9Q{Xj)) is 0-Hermitian. The following Lemma summarizes the consequences. 



Lemma 2.3 Let A be normal and 9 such that (11) is satisfied. Then there is a it £ 



such that Q {A) = i ■ 9tt{9Q{A)), and thus the 9-Toeplitz decomposition of A can be written 
as 

A = e{A) + i-97r(9e{A)), (14) 

or, in terms of N = 9 A and its classical Toeplitz decomposition, N = %{N) + i ■ tt{'H{N)). 
The eigenvalues of A lie on the image C of the function 

c:R^C: p ^ 9p + i ■ 9tt{p). (15) 

Proof. Apply the spectral theorem for normal and Hermitian matrices. □ 

Remark 2.4 Note that A is essentially Hermitian if and only if there exists a G T such 
that the interpolating polynomial vr G 7^^(M). In fact, vr may then even be in V^{I). 

Remark 2.5 If for some normal matrix A the degree of the interpolation polynomial equals 
two for some value of 9, it may well be of degree n — 1 for almost all other values of 9, since 
then the eigenvalues lie on a rotated parabola, as depicted in Figure [T| There does not seem to 
be an easy way to determine for which the polynomial degree is minimal, although for each 
given 6 the polynomial can be computed exactly in a finite number of arithmetic operations 
without knowing the eigenvalues. 



2.2.2 Computing the polynomial vr for given 6 & T 

For almost any fixed value of 9, the interpolation polynomial vr belonging to a normal matrix 
A can be computed, without knowing the eigenvalues of A, in a finite number of arithmetic 
operations, providing us with a curve C C C on which all eigenvalues of A lie. Indeed, if the 
degree of vr equals k then A — Q{A) is a linear combination of 

i9I, i9Q{A), i{9e{A)f, i{9Q{A))^. (16) 
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Figure 1: Interpolating polynomials of different degree for different values of 9. The seven 
asterisks represent the eigenvalues. 

Making the combination explicit is equivalent to finding the coefficients of tt. To obtain vr in 
practice, notice that for any v G C", 

Tr(9Q{A))v G IC''(eQ{A),v) = span{v, 9e{A)v, (ee{A))''v}, (17) 

the Krylov subspace for 6Q{A) and v. Since 9Q{A) is Hermitian, an orthonormal basis for 
}C''{6Q{A),v) can be constructed using a three-term recursion, and solving the linear system 
can be done cheaply. These, and other considerations, led Huhtanen to the development 
of efficient structure preserving eigensolvers [8] and linear system solvers [9] for problems 
involving normal matrices. 

2.3 Commuting normal matrices 

Because we will need to draw conclusions about the eigenvalues and eigenvectors of commuting 
normal matrices, here we recall a well-known result and give its complete proof. 

Lemma 2.6 Normal matrices A and E commute if and only if they are simultaneously uni- 
tarily diagonalizable. 

Proof. If both A = W*AW and A = W*EW are diagonal for a unitary matrix W, then 
clearly [A,E] = M^[A,A]VF* = 0. Conversely, assume that [A,E] = 0. Let [/ be a unitary 
matrix such that A = U*EU is diagonal with multiple eigenvalues being neighbors on the 
diagonal of A. Thus 

A = 
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where rrij denotes the multiphcity of 6j. Let S = U*AU and write Spq for its entries. Then 
equating the entries of SA and AS in view of the relation 

[S, A] = U*[A,E]U = 0, 

shows that Spq = whenever dp ^ 5q. Thus, S is block-diagonal with respective blocks 
Si, . . . , Si of sizes mi, . . . , mi. For each j £ {!,...,£}, Sj is normal. Let Qj be unitary and 
such that Aj = Q*-SjQj is diagonal. Then with 



Q 



Qi 



Qi 



and A 



Ai 



and W = UQ 



(18) 



we find that A = Q*AQ = W*EW and A = Q*SQ = W*AW. 



□ 



Remark 2.7 In case all eigenvalues of E are distinct, then S itself is already diagonal and 
the proof is finished. In case E has an eigenvalue, say 6i, of multiplicity mi > 1, then there is 
freedom in the choice for the first mi columns ui, . . . , Umi of U that correspond to 6i. Even 
though each choice diagonalizes E, not each choice diagonalizes A as well. This is expressed as 
S having a diagonal block 5*1 of size Till. Writing for the matrix with, coluirms tii, • ■ ■ ; 
we have that 

AUi = UiSi, (19) 

hence the column span of Ui is an invariant subspace V of A. The matrix Qi determines, 
through the transformation Wi = UiQi, an orthonormal basis for V of eigenvectors of A. If 
^i has multiple eigenvalues, again there may be much freedom in the choice for Qi. 



3 Normality preserving augmentation 



In this section we revisit, from an alternative point of view, a problem studied by Ikramov 
and Eisner in [11]. It concerns the augmentation by a number m of rows and columns of a 
normal matrix in such a way, that normality is preserved. Our analysis differs from the one 
in [11], and we add details on the eigendata of A^^- in terms of those of A. 



Normality preserving augmentation. Let A £ M"^"(C) be normal. Characterize all 
nxm matrices V,W e M"^™(C), and all Y G M™^™(C) such that A+, where 



A, 



A 

W* 



V 
Y 



(20) 



is normal, too. In other words, characterize all normality preserving augmentations of A. 



Remark 3.1 Note that Hermicity, 0-Hermicity, and essentially Hermicity preserving aug- 
mentation problems are all trivial, because each of these properties is inherited by principle 
submatrices. For unitary matrices, this does not hold. However, if A and are both unitary, 
their rows and columns all have length one. Thus V = W = and also Y is unitary. This 
trivially solves the unitarity preserving augmentation problem. When we consider normality 
preserving augmentation, matters become less trivial. 
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3.1 The normality preserving augmentation problem for m = 1 



First consider the case m = 1, and write v^w^y instead of V,W,Y. It is easily verified that 
= if and only if 



WW 



vv 



w*w = v*v, and A*v + yw = Aw + yv. 



(21) 



The two leftmost relations hold if and only if w = (pw for some </> G T. The rightmost relation 
may add further restrictions on v,w and (f). Before studying these, note however that (j) is the 
square of a unique 9 G Tu C T, where 



Tu = {tGT I arg(r) e [0,7r)}. (22) 
')w that better reveals the underlying structure, which is 
bw 44> V = f'w 44> 9v = 6w = u V = 9u and w = 6u (23) 



This yields a reformulation of v 

32 



for some u G C". Further restrictions on u and 6 G Tu follow from substituting v = 9u and 
w = Ou into the rightmost equation in (21), which after some rearrangements results in the 
requirement that 

e^{A)u = e^{y)u. (24) 

we have solved the 



2.2 



Because eigenpairs of Q (A) were already characterized in Corollary 
augmentation problem for m = 1 . The below theorem summarizes this solution constructively 
in terms of the eigendata of A. Note that this result was proved already in [llj, though in a 
different manner. 

Theorem 3.2 Let y £ C, and £ : M ^ C : p y + p9 a line in C through y with slope 
9 £ T\j. Let Ai, . . . Ap be the eigenvalues of A that lie on i. Then, the matrix 



A, 



A V 
w* y 



(25) 



is normal if and only if v = 9u and w = 9u, where u is a linear combination of eigenvectors 
corresponding to Ai, . . . A„, with the convention that n = if p = 0. 



Proof. Corollary 



2.2 



shows the relation between eigendata of Q (A) and A, and together 
with the derivation in this Section this proves the statement. □ 



3.2 The eigenvalues of the augmented matrix 

Next, we augment the analysis in |11| with a study of the eigenvalues of A^ in relation to 
those of A. Let Ap G M^^^(C) be the diagonal matrix whose eigenvalues are the p eigenvalues 
of A that lie on ^ : M — )• C : p I— 7- y + p9. Then, as already mentioned in the proof of Corollary 



2.2 



A. 



.p-9R + yLp (26) 

for some real diagonal matrix R. Let U G M"^"(C) be any unitary matrix whose last p 
columns are eigenvectors of A belonging to Ai, . . . , Ap and let Up contain those last p columns 
of U. Then, assuming that A and A^ are normal. Theorem 



3.2 



shows, with r 



UpU, that 



" u 


* 


A 


V 




' u 


1 




w* 


y . 




1 



B 



A, 



9r 

y 



(27) 
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Moreover, 



9r* 



Or 

y 



9Rj^ + yIp+i, where Rj, 



R 

r* 



(28) 



The above observations reveal some additional features of the solution of the augmentation 
problem, that we will formulate as another theorem. 

Theorem 3.3 The only normality preserving 1- augmentations of A are the ones that, on an 
orthonormal basis of eigenvectors of A, augment a p x p essentially Hermitian submatrix of 
A. Hence, n — p eigenvalues of A are also eigenvalues ofA^. The p + 1 eigenvalues of Aj^ 
that remain lie on the same line as, and are interlaced by, the remaining p eigenvalues of A. 



Proof. The block form in (27) shows that the eigenvalues of B are eigenvalues of both A 
and Aj^ , whereas ( 28 ) shows that to locate the remaining p + 1 eigenvalue of A. 
needs to observe that the real eigenvalues of R interlace those of 



one only 
□ 



Remark 3.4 Note that the case p = 0, covered by Theorem 3.2 , is also included in the above 
analysis if one is willing to interpret on the same line as the remaining p eigenvalues of A as 
any line in C. This just reflects that the additional eigenvalue y S C of A^ can lie anywhere. 



For an illustration of the constructions of Theorems 3.2 and 3.3 see Section 6.1 There, 



we will augment a given 3x3 matrix in two different ways and compute the eigenvalues of 
the augmented matrix. 



3.3 Normal matrices with normal principal submatrices 

By applying the procedure for m = 1 several times consecutively, also m-augmentations with 
m > 1 can be constructed. In particular, all the normal matrices having the property that all 
their principal submatrices are normal, can be constructed. Since generally, normal matrices 
do not have normal principal submatrices, this shows that the m-augmentation for m > 1 
has not yet been completely solved. In Section [5] we will study the principal submatrices 
of normal matrices from the point of view of Section 2.2, This will also give more insight 
into the m-augmentation problem for m = 2, that already proved to be very complicated in 
In particular, we will give a procedure to augment A that does not reduce to m-fold 
application of the 1-augmentation. Before that, we investigate normality preserving normal 
perturbations. Apart from being of interest on its own, it will be needed in Section [5] 



4 Normality preserving normal perturbations 

In this section we consider a question related to the augmentation problem, and we will study 
it using the same techniques. In particular, we will study normality preserving ^-Hermitian 



perturbations, that play a role in the augmentation problem in Section 5.1 



Normality preserving normal perturbation. Let A € M"^"(C) be normal. Character- 
ize all normal E such that = A-\-E is normal. In other words, characterize the normality 
preserving normal perturbations E of A. 

Since A = Q{A) + Q-^{A), and both 0(^4) and Q^{A) are normal, this shows that the sum 
of normal matrices can be literally any matrix and thus that the above problem is non-trivial. 
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4.1 Generalities 

To start, we formulate a multi-functional lemma that summarizes the technicalities of writing 
out commutators of linear combinations of matrices. 

Lemma 4.1 Zei yl, £^ G M"^"(C) and^,^^C Then with 9 = 

[jA + fiE,{-fA + i^Ey] = \jf[A,A*] + 2j]I-e{[A,E*]) + \fi\^[E,E*], (29) 

and 

e{[A, E*]) = [e{A),e^{E*)] + [e^{A), e{E*)]. (30) 

Therefore, if A and E are normal, then jA + fiE is normal if and only if 

ei[A,E*]) = 0, (31) 
or in other words, if and only if [A, E*] is ^fi-skew-Hermitian. 

Proof. Straight-forward manipulations with the commutator give the statements. □ 

Corollary 4.2 Let A, E £ M"^"(C) be normal. Then A + E is normal if and only if jA + 
IjlE + al is normal for alJ 7, /u, a G C with the restriction that 'jjl G M. 

Corollary 4.3 IfA,Ee M"^"(C) are normal then [A, E*] = ^ [A, E] = 0. Thus if either 
term vanishes, both A + E* and A + E are normal. 



Proof. As was shown in the proof of Lemma 2.1, E and E* have the same eigenvectors. 



Thus, A and E* are simultaneously unitarily diagonalizable if and only if A and E are. Lemma 



2.6 now proves that [A, E*] = 44> [^4, E] = 0, and Lemma 4.1 proves the conclusion. □ 



Corollary 4.4 The sum of ^A and jiE with £ C and A and E Hermitian is normal if 
and only if^n GM or [A, E] = 0. 

Proof. The commutator of Hermitian matrices is always skew-Hermitian. Thus for [^4, E] 



to be 0-Hermitian in (31), jj-i must be real, or [j4,i?] should vanish. □ 



4.2 Normality preserving normal rank one perturbations 

This section aims to show the similarities between the normality preserving normal rank one 



perturbation problem and the m-augmentation problem of Section 3.1 Indeed, let E = vw* 
with v,w G C". Then E is normal if and only if 

llwll^fw* = llvlpiiitD*, (32) 

and thus if and only if ?; = zw for some z G C Write z = 9p with ^ G T and p G M, p > 0. 
This shows that a rank-1 matrix E is normal if and only if E is ^-Hermitian, 

E = 9uu*, 9eT. (33) 

With A G M"'^"(C) normal, we will look for the conditions on u G C" and G T such that 



A + 9uu* is normal. Since 9uu* is 0-Hermitian, by (30) in Lemma 4.1 



[{A + 9uu*),{A + 9uu*)*] = 29[Q^{A),uu*], (34) 
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thus, A + duu* is normal if and only if Q-^{A) and uu* commute. According to Lemma 



2.6[ this is if true and only if they are simultaneously unitarily diagonalizable. For this, it 

we will 



3.2 



is necessary and sufficient that u is an eigenvector of G (^). As in Theorem 
formulate this result constructively in terms of the eigendata of A. 

Theorem 4.5 Let y £ C, and i : M ^ C : p y + p6 a line in C through y with slope ^ G T. 
Let Ai, . . . Ap be the eigenvalues of A that lie on i. Then, the matrix 



A+ = A + Ouu* 



(35) 



is normal if and only if u is a linear combination of eigenvectors corresponding to Ai, . . . Ap, 
with the convention that u = if p = 0. 



Proof. Corollary 2.2 shows the relation between eigendata of Q (A) and A, and together 
with the derivation above this proves the statement. □ 

Remark 4.6 An interesting consequence of adding the normality preserving normal rank 
one perturbation E = Ouu* is that 



e^{A + 6uu*) = e^{A) + e^ 



1UU 



(36) 



because Ouu* is 0-Hermitian. Thus, the conditions under which adding another 0-Hermitian 
rank one perturbation F = Oww* to A + E lead to a normal A + E + F are identical to the 



conditions just described for E. We will get back to this observation in Section 4.4 
4.3 The eigenvalues of the perturbed matrix 

To study the eigenvalues of A^ in relation to those of A, let Ap G Mp^p{C) be the diagonal 
matrix whose eigenvalues are the p eigenvalues of A that lie on i : M ^ C : p y + pO. Then 



Ap = OR + yLp 



(37) 



for some real diagonal matrix R. Let U G M"'^"(C) be any unitary matrix whose first p 
columns are eigenvectors of A belonging to Ai, . . . , Ap. Then, assuming that A and A^ are 



normal. Theorem 4.5 shows that 



U*A+U = U*{A + Ouu*)U 







Orr* 




+ 


B 





(38) 



because n is a linear combination of the first p columns of U . This leads to the following 
theorem, in which we summarize the above analysis. 

Theorem 4.7 The only normality preserving normal rank one perturbations of A are the ones 
that, on an orthonormal basis of eigenvectors of A, are 0-Hermitian rank one perturbations 
of a p X p 0-Hermitian submatrix. Hence, n — p eigenvalues of A are also eigenvalues of A^. 
The remaining p eigenvalues of A^ are the eigenvalues of 

Ap + Orr* =0{R + rr*) + yLp, (39) 

and these interlace the p eigenvalues of A on i with the additional p + 1-st point +oo0. 
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Proof. The eigenvalues of are the eigenvalues of B together with the eigenvalues of the 
matrix in (39). Obviously, all eigenvalues of B are eigenvalues of A as well. Since rr* is a 
positive semi-definite rank one perturbation of R, the eigenvalues pi < ■ ■ ■ < Pp oi R + rr* 
and the eigenvalues ri < ■ ■ ■ rp of R satisfy 

< yOi < r2 < • • • < pp-i <rp< pp (40) 

as a result of Weyl's Theorem [12] . Multiplying by 9 and shifting over y yields the proof. □ 



Remark 4.8 Note that if p = 1, only one eigenvalue is perturbed, and pi = ri + In 
terms of the original perturbation E = Ouu* this becomes A = A + where A is the 

eigenvalue of A belonging to the eigenvector u. 

As a consequence of the following theorem, it is possible to indicate where the eigenvalues 
of the family of matrices A + tE are located. This can only be done for normal normality 
preserving perturbations. 

Theorem 4.9 Let A,B e M"^"(C) be normal. Consider the line £ through A and B, 

^:M^M"^"(C): t^tA+{l-t)B. (41) 

If E = B — A is normal, all matrices on i are normal; if E is not normal, A and B are the 
only normal matrices on I. 



Proof. Observe that £(t) = A + (1 — t)E. If E is normal, then Corollary 4.2 shows that all 
matrices ^A + pE with 771 S M are normal, which includes the line I. Assume now that E is 
not normal. Because A and B = A + E both are normal, (29) in Lemma 4.1 gives that 

{l-t){2'H{[A,E*]) + {l-t)[E,E*]) = with [E,E*]^0. (42) 

The solution t = 1 confirms the normality of A, and the linear matrix equation 

2n{[A,E*]) + {l-t)[E,E*]=0 with [E,E*]^0 (43) 

allows at most one solution in t which, by assumption, is t = 0. □ 

Thus, any line in M"'^"(C) parametrized by a real variable that does not lie entirely in the 
set of normal matrices, contains at most two normal matrices. 



Remark 4.10 Lemma 2.1 shows that if E £ M"^"(C) is such that A + E is normal, then 

a{A + E)C a{Q{A) + G(^)) x a{Q^{A) + e^(^)), (44) 

and perturbation theory for 0-Hermitian matrices can be used to derive statements about the 
eigenvalues oi A + E. According to Theorem 4.9, if E itself is normal too, this relation is 
valid continuously in t along the line A + tE: 

a{A + tE) C a{e{A) + te{E)) x a{Q^{A) + tQ^{E)). (45) 



For non-normal E this would generally not be true, as will be illustrated in Section 6.2.1[ 
Corollary 4.11 As a result of Theorem \4.S\ and Remark \4.10[ the perturbed eigenvalues of 

t^A + tOuu*, 0<t<l, (46) 
seen as functions of t, are line segments that all lie on the same line with slope 6. 
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4.4 ^-Hermitian rank-A; perturbations of normal matrices 



In Section 5.1 we will encounter normality preserving 0-skew-Hermitian perturbations. In 
this section we will fully characterize those. For this, consider for given k with 1 <k <n the 
0-Hermitian rank-Zc matrix 

E = OH, where 6* E Tu and H = H* . (47) 



Let A be normal. Then, since Q (E*) = 0, Lemma 4.1 shows that A + E is normal if and 
only if 

[e^{A),H] = 0. (48) 



By Lemma 2.6 this is equivalent to @ (A) and H being simultaneously diagonalizable by a 
unitary transformation U. Thus, H needs to be of the form 

H = UAU* (49) 

where A E M^^'^(]R) is diagonal with diagonal entries 6i, . . . ,Sk and the columns ui, . . . ,Uk 
of U are orthonormal eigenvectors of @-^{A). But then, writing 

E = Ei-\ \- Ek, where for ah j E {1, . . . , k}, Ej = OdjUjU*, (50) 



the observation in Remark 4.6 reveals that perturbing ^ by is equivalent to perturbing A 
consecutively by the rank one matrices Ei, . . . , E^. The above analysis is summarized in the 



following theorem. An illustration of this theorem is provided in Section 6.2.2 



Theorem 4.12 Let E = OH he a 6-Hermitian rank k perturbation of a normal matrix A. 
Then E is normality preserving if and only if E can be decomposed as 

E = Ei + --- + Ek, (51) 

where Ei, . . . , E^ are all normality preserving 9-Hermitian rank one perturbations of A. In 
fact, for each permutation a o/ {1, . . . , /;;} and each m E {1, . . . , A;}, the partial sum 

m 

A + Y.^'^U) (52) 
i=i 

is then normal, too. 



Remark 4.13 In accordance with Remark 2.7, if A in (49) has multiple eigenvalues, there 
exist non-diagonal unitary matrices Q such that A = QAQ*. As a result, H can be written 
as ZAZ* where the orthonormal columns of Z span an invariant subspace of Q'^{A). This 
implicitly writes the perturbation 9H as a sum of rank-1 normal perturbations that do not 



necessarily preserve normality. This aspect is also illustrated in Section 6.2.2 



Theorem 4.14 The only normality preserving 6-Hermitian rank-k perturbations of A are the 
ones that, on an orthonormal basis of eigenvectors of A, are 9-Hermitian perturbations of 9- 
Hermitian submatrices of size si x si, . . . , x Sm of ranks ki, . . . , km, where ki-\-- ■ ■ + km = k. 
As a result of this perturbation, at most si + • • • + eigenvalues of A are perturbed, which 
are located on at most m distinct parallel lines ii, . . . , im, defined by 9 £T and yi, . . . , E C 
as 

ij-.^^C: p^yj + 9 p. (53) 

Moreover, the eigenvalues of A -\- tE, t E [0, 1] connect the perturbed eigenvalues of A with 
those of A by line segments that lie on £i, . . . , Im- 
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Proof. Write 9H = 9{5iuiu\ + • • • + 6kUku\), where ui, . . . ,Uk are eigenvectors of Q-^{A), 



and repeatedly apply Theorem 4.7 The statement about the eigenvalues of A + tE follows 
from Corollary 4.11 □ 



For a qualitative illustration of the effect on the eigenvalues due to a rank-1 0-Hermitian 
perturbation and a rank- A; 0-Hermitian perturbation of a normal matrix, see Figure [2| The 
asterisks in the pictures are eigenvalues of A, the circles represent different choices for y, and 
the boxes are the perturbed eigenvalues. 




(3 S(2) 



3t{z) 




Figure 2: Eigenvalue perturbation by a rank-1 matrix (left) and a rank-/c matrix (right). 



Remark 4.15 If H in (47) is semi-definite, then the eigenvalues oi A + tE all move in the 



same direction over the parallel lines £i 



from Theorem 4.14 



4.5 Normality preserving normal perturbations of normal matrices 

The above analysis of 0-Hermitian normality preserving perturbations also gives sufficient 
conditions for when normal perturbations E of the form E = OiHi + O2H2 with Hi,H2 
Hermitian and 01,02 £ T are normality preserving in case 02 does not equal it^i . Notice that 
in order for E to be normal itself, [Hi,H2] = by Corollary 4.4 Obviously, E is in general 
not 0- Hermitian for some value of 0. Nevertheless, the following holds. 

Theorem 4.16 Let A be normal, Hi,H2 Hermitian with [Hi,H2] = 0, and 0i,02 S T with 
0102 Then E = Ei + E2 is a normality preserving normal perturbation of A if Ei and 
E2 both are normality preserving perturbations of A. 

Proof. Corollary |4.4| covers the normality of E. Furthermore, assuming that A + Ei is 
normal, (48) gives that A + Ei + E2 is normal if and only if [9-'-(j4 + £'1), H2] = 0, where 0-*- 



indicates the 02-skew Hermitian part. But 

[e^(^ + Ei),H2] = [Q^{A),H2] + [e' 
because [Hi,H2] = 0, proving the statement. 



hHi),H2] = [e^{A),H2], 



(54) 
□ 



Remark 4.17 A similar result holds for normal perturbations E = Ei + ■ ■ ■ -\- E^ where 
Ej = 9jHj with Hj Hermitian and 9j G T for each j G {1, . . . , A;}. Moreover, each normal 
perturbation E can be written in this form in several different ways. 

We are now ready to return to the augmentation problem of Section [4] and to study 
augmentations of normal matrices A that result in an augmented matrix A^ whose eigenvalues 
all lie on the graph of a quadratic polynomial. 
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5 Further augmentations 



We now return to the m-augmentation problem of Section |4j We concentrate on the case 
m > 1 and on augmentations yl+ that have some additional structure in the spirit of Section 



2.2 For instance, if all eigenvalues of lie on a line, then is essentially Hermitian, a 
property that is inherited by principal submatrices. In that case it is clear which matrices 
A can be augmented into A^. The next simplest case, if simple at all, is the case where all 
eigenvalues of A^ lie on the image of a quadratic function in a rotated complex plane. 



5.1 Augmentations with all eigenvalues on a quadratic curve 

Assume that for G T there exists a quadratic vr G T'^(M) such that 

A+ = e{A+) + i-9TT {eG{A+)) , where A+ - 



A 

W* 



V 
Y 



(55) 



Then A^ is normal with all its eigenvalues on the rotated parabola C C C defined as the 
image of g, where 

q:M^C: p^9p + i- 9tt{p) (56) 

and 



tt{x) = ro + rix + r2X with rQ,ri,r2 G 



(57) 



Remark 5.1 Throughout this section we assume, without loss of generality, that r2 > 0. 
The case r2 = 0, as argued above, is trivial and concerns essentially Hermitian matrices, 
whereas the case r2 < can be avoided by replacing 9 by —9, which is nothing else than a 
trivial change of coordinates that transforms the polynomial vr into — vr. 



The curve C now divides the complex plane C in three disjoint parts 

c = c+ucuc_, 

where C+ is the open part of C that lies on the one side of C that is convex. 

5.1.1 The principal submatrices of and their eigenvalues 

For convenience, write 

X = 9e{A), M = 9e{Y) and 2Z = 9V + 9W, 

then 



(58) 



9e{A^ 



X 

z* 



z 

M 



and {9QiA+)f 



X^ + ZZ* XZ + ZM 
Z*X + MZ* M'^ + Z*Z 



(59) 



(60) 



with the block form of A^ displayed in ( 55 ) yields that 



Thus, explicitly evaluating vr at A+ using the block forms in (|60|), and comparing the result 

(61) 



A = e{A) + i ■ 9tt{9Q{A)) + i ■ 9r2ZZ* 



and 



Y = Q{Y) + i ■ 97r{9e{Y)) + i ■ 9r2Z*Z. (62) 
Results that follow will sometimes be stated for A only, even though similar statements 



obviously hold for Y . The first proposition simply translates (61) in words. 
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Proposition 5.2 The n x n leading principal suhmatrix A of Aj^ is a 9-skew-Hermitian 
rank-k perturbation of a normal matrix that has all its eigenvalues on C, and k < min(m,n). 



Lemma 5.3 If A in (61) is normal, then (j{A) C C UC+. 



Proof. If A is normal, then obviously i ■ 9r2ZZ* is a normality preserving ^-skew-Hermitian 
perturbation of the normal matrix Q{A) + i ■ 67r{9Q{A)) that has all its eigenvalues on C. By 



Theorem 4.14, each eigenvalue A G C that is perturbed, moves along the line £ : R - 
/9 I— )• A + i • 9p. Note that i is vertical in the 0-rotated complex plane. By Remark 4.15 since 



ZZ* is positive semi-definite, the direction is the same for each perturbed eigenvalue, and is 



determined by the sign of r2. In Remark 5.1 we assumed that r2 > 0, and thus the direction 



is directed into C-(- defined in (58). □ 



Remark 5.4 Note that a multiple eigenvalue A of Q{A) + i ■ 9tt(9Q{A)), located on C, may 
be perturbed by i ■ 9r2ZZ* into several distinct eigenvalues of A. Those will all be located on 
i:R^C: X + i-9p with p > 0. 



Corollary 5.5 Assume that A in (61) is normal. Then cr{A) C C if and only is block 
diagonal with blocks A and Y . 



Proof. If Z ^ then tTace{ZZ*) ^ and at least one eigenvalue is perturbed. Lemma 5.3 



shows that a perturbed eigenvalue cannot stay on C and necessarily moves from C into C+. □ 
5.1.2 Augmentations with eigenvalues on a quadratic curve 



The observations in Section 5.1.1|can be reversed in the following sense. Given j4, we choose 



a parabolic curve C and construct Z G M"^™(C) such that i-9r2ZZ* perturbs the eigenvalues 
of A onto C. Then we use C to define the corresponding m-augmentation of A. 

Corollary 5.6 Necessary for a normal matrix A to be m-augmentable into a normal matrix 
A^ with all eigenvalues on a quadratic curve C is that cr(A) C C UC+. 



Proof. This is just another corollary of Lemma 5.3, □ 



Clearly, for any given finite set of points in C, there are infinitely many candidates for 
such quadratic curves C. It is the purpose of this section to show that each of this candidates 
can be used, and to construct essentially all possible corresponding augmentations A^. 

Theorem 5.7 Let A G M"^"(C) be normal, and let 9 G T and vr G P^^M) be such that 
cr{A) C C U C_|_, where C is the graph of 

q:R^C: p ^ 9p + i ■ 9tt{p). (63) 

Then there exist p- augmentations A^ of A such that 

a{A+) C C, (64) 

where p is the number of eigenvalues of A in Cj^. 
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Proof. Write Ap G MP^^(C) for the diagonal matrix with precisely the eigenvalues Ai, . . . , Ap 
of A that do not lie on C and let Up € M"^p(C) have corresponding orthonormal eigenvectors 
ui, . . . ,Up as columns. Since cr{Ap) C C+, for each j £ {1, . . . ,p} there exists a positive real 
number such that 

Xj-i- Bij G C. (65) 
Write H G MP^p(C) for the diagonal matrix with cr(Hp) = • • ■ , \/%>} and set 



Z — Up^p. 



(66) 



By Lemma 2.1, the columns of Up are also eigenvectors of Q{A) and thus i ■ 6ZZ* is a 



0-skew-Hermitian normality preserving perturbation of A. Moreover, 

u*{A-i- ezz*)uj = x-i-ecec. 



(67) 
QiA) 



Because A — i ■ 9ZZ* is normal with all eigenvalues C, the equality Q{A — i ■ 6ZZ*) 
leads to 

A-i-0ZZ* = 9@{A) +i-dTT(eQ{A)). (68) 



(69) 



The assumption r2 > 0, justified in Remark |5.1[ now gives that 

A = 9e{A) + i ■ eTT(eQ{A)) + i-er2 



and according to (61 ) this is precisely the nxn leading principal submatrix of A+, where A 

ZQ 



is defined as A^ = 9H + i ■ 9tt{H), where 

H = 



9e{A) z 
Z* M 



with Z 



and M,Q e 



*(C) are arbitrary Hermitian and unitary matrices. 



(70) 



□ 



For a given nxn normal matrix A, the typical situation is that after selecting 9 suitably, 
at least p of its eigenvalues do not lie on a quadratic curve, and a matrix Z of rank p is needed 
to push those outliers onto C. In Figure 3, already 3 of the 7 eigenvalues of A, indicated by 
asterisks, lie on the quadratic curve C, and a rank-4 matrix Z is needed to push the remaining 
4 eigenvalues onto C. After that, the augmented matrix can be formed. 

Remark 5.8 It is, of course, possible to move each eigenvalue of A from C_|_ onto C as a result 
of an arbitrary amount of rank-1 perturbations. This would increase the number of columns 
of Z, and give m-augmentations of A with m > p. However, this would not increase the rank 



of Z, and ZZ* would remain the same. Together with the analysis of Section 4.4 that shows 



which 0-Hermitian perturbations are normality preserving, this shows that in essence, each 



p-augmentation A^ with Z of full rank, is of the form (70). In Section 6.3 we give an explicit 



example of the construction in the proof of Theorem 5.7 
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Figure 3: Illustration of the construction in the proof of Theorem 5.7, Already 3 of the 7 
eigenvalues of A lie on the quadratic curve C, and a rank-4 matrix Z is needed to push the 
remaining 4 eigenvalues onto C. 



5.1.3 Augmentations without computing eigenvalues 

So far, explicit knowledge about the eigenvalues and eigenvectors of A was used to construct 
augmentations A-^-. There are, however, cases in which it is sufficient to know the polynomial 
curve £ on which the eigenvalues of A lie. To see this, assume that A is normal and 

A = e{A)+i-e7r(eQ{A)), 7rer^''{R) and tt ^ V^^-^R) (71) 

for some integer k > 1. Since vr has even degree, there exist polynomials p G 'P^(]R) such that 

7r(x) - p{x) > for aU x G M. (72) 

This implies that the matrix (vr — p){9Q{A)) is positive semi-definite, and hence it can be 
factorized as 

{7r-p)(0e{A)) = ZZ*, (73) 

after which we have that 

A = Q{A)+i-9p(9e{A)) + i-9ZZ*. (74) 

It is trivial that i ■ 6ZZ* is a normality preserving perturbation of Q{A) + i ■ 6p{9Q{A)), and 
by choosing between 9 and —9, as explained in Remark 5.1, this leads to an m-augmentation 
of A, with generally m = n — 1. Section [2.2.2| explained how vr can be computed in a finite 
number of arithmetic operations, and the same is valid for the factorization (73). Of course, 
the problem of finding a minorizing polynomial p G 'P(M) may prove to be difficult in specific 
situations. 
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5.2 Polynomial curves of higher degree 



If one tries to generalize the approach of Section 5.1 to polynomial curves C of higher degree, 
the situation becomes rapidly more complex. As an illustration, consider the cubic case. The 
third power of the matrix 9Q{A^) in (60) equals 



+ XZZ* + ZZ*X + ZMZ* X'^Z + XZM + ZM"^ + ZZ*Z 
Z*X'^ + Z*ZZ + MZ*X + M'^Z + Z*ZM + MZ*Z + Z*XZ 



(75) 



and thus, comparing the leading principal submatrices, 

A = Q{A) + i ■ eT:(eA) + i ■ Br^ZZ* + i ■ er^{XZZ* + ZZ*X + ZMZ*), (76) 

where ra is the coefficient of of vr. Thus, ^ is a rank-/c, with k < 2m, ^-skew-Hermitian 
perturbation of the normal matrix 0(^4) + i • 9tt{6A). Of course, if ZZ* commutes with 0(^), 
then it commutes with X, and this may help the analysis. However, it becomes much harder 
to control the perturbation in such a way, that A will be augmented into a matrix j4+ with 
a{A^) C C. Therefore, we will not pursue this idea in this paper. 



6 Illustrations 

In this section we present some illustrations of the main constructions and theorems of this 
paper. By making them explicit, we hope to create more insight in their structure. 

6.1 Illustrations belonging to Section [3] 



This example illustrates Theorems 3.2 and|3.3[ Let A be the matrix 



A 



2i 
2 + i 




(77) 



Take y = 1 and choose ii to be the line through y = 1 and the eigenvalue 2 + i of A, 

l + i 



'1 : 



pi-^ l + 9p, with 6 



V2 



(78) 



Thus, in order for A^ to be normal, u must be a multiple of 62, and v = Ou, w* = 6u* , which 
yields that for all /.i G C, 



(79) 





' 2i 
















2 + i 













-3 










eji 





1 



is a normal augmentation of A. Moreover, for y = 1 and this choice for ii, these are all the 
normal augmentations of A. The eigenvalues of A+ that are not eigenvalues of A are the 
eigenvalues of 



2 + i Op ' 




" 1 + y/2e 0p 


= d 


' V2 p' 


Op 1 




eji 1 




p 



+ 1, 



(80) 



19 



and thus equal to 

x = e(^^±V2fi^ + ij +1, (81) 

which He on ii and have the eigenvalue 2 + i that was perturbed as average, as is depicted in 
the left picture of Figure |4j The stars in both pictures are eigenvalues of A, the open circles 
represent different choices for y. Perturbed eigenvalues are indicated by squares. 




Figure 4: Left: the choice y = 1 and the line li trough the eigenvalue 2 + i. Right: the two 
lines £2 and £3 are depicted that correspond to the choice y = —2 + 3i. 

A second option is to choose y = — 2 + 3i instead oi y = 1. This gives other possibilities 
to construct normality preserving augmentations. The first one is to choose £2 through y and 
—3, which shows that u must be a multiple of 63 and v = 9u, w* = 9u where 6 = (l + 3i)/-v/l0, 
which is a similar situation as for y = 1. The second non-trivial option is to choose £3 through 
y and both 2i and 2 + i. Then with 6 = (— 2 + i)/\/5, we may take u as any linear combination 
of ei and 62, showing that 





' 2i 








6*7 


A+ = 





2 + i 





Ofi 








-3 







_ Oj 


eji 





-2 + 3i 



(82) 



is normal for all 7, /x G C. For the given value of y those two options are the only possible 
ways to construct normality preserving augmentations. See the right picture in Figure |4] 



6.2 Illustrations belonging to Section |4] 

6.2.1 Normality preserving non- normal perturbation 



First we illustrate Theorem 4.9 by presenting an example of a line through two normal matrices 
that contains only two normal matrices. For this, let 



1 

1 



B 



1 
-1 



and E = B - A 




-2 



Thus, ii^ is a non-normal normality preserving perturbation of A. Hence, apart from A and 
B, no other matrix of the form A + tE with t € M is normal. Indeed, 



[A + tE, {A + tE)* 



1 - (1 -2t)2 






(1 - 2tf - 1 
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and this matrix is only zero for t = and t = 1. Moreover, the eigenvalues of A + tE are 

Vl - 2t and - Vl^^, 



whereas the corresponding sums of the eigenvalues of ^{{A + tE) and 5(^4 + tE) equal 

(1 - t) + it and - (1 - t) - it. 

Thus, for instance at t = ^, the eigenvalue zero of A + tE is not the sum of the eigenvalues 
of the the Hermitian and skew-Hermitian parts of A + tE. 



6.2.2 Normality preserving 0-Hermitian perturbations 



Next, we illustrate Theorem 4.12 This concerns normality preserving 6 -Hermitian perturba- 
tions. For this, we take 9 = 1 and consider the matrix A-\- E, where 



A 



1 
i 
1 + i 



and E = uiu* + 2u2U2 



with ui,U2 S mutually orthonormal. Then A + E is normal if and only if both ui and U2 
are eigenvectors of 

[000 

2e^{A) = 2S{A) = i 

i 

However, A-\- E, where E = uiul + U2U2 with ui,U2 G C'^ mutually orthonormal, is normal 
if and only if both ui and U2 are linear combinations of the same two eigenvectors vi and f 2 
of 2S{A). Thus, with 



Ul 



V2 
1 
1 



U2 



-V2 
1 
1 



where vi 



1 





and V2 



and El = uiu^ and E2 = U2U2 we have that 

2 V2 V2 
\/2 1 1 + 
V2 1 1 



El + E2 



1 

16 



1 

16 



2 

-V2 
-V2 



-V2 
1 
1 



'V2 
1 
1 



1 

16 



4 







IV2 

kV2 




2 2 
2 2 



is a normality preserving rank two perturbation of A, written as the sum of two rank one 
normal perturbations Ei and E2 that each does not preserve normality. However, we also 
have that 



-^VlVi + -V2V2 = - 



and this expresses the same perturbation E as a sum of two normality preserving normal 
rank-1 perturbations. The eigenvalues of A -\- tE are 1 + |t due to the term jVif*, together 
with the eigenvalues of 



" 4 





" 


1 

+ 16 


" 





" 


1 

~ 16 


" 4 





" 


















2 
2 


2 
2 






2 
2 


2 
2 



' i ■ 


1 


" 1 1 " 


1 + i 


+ -t 
8 


1 1 



which are i + 



1 



1 + 4 + 



1 + 



16 



(83) 
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which are due to the term \v2V2 



As stated in Theorem 



4.14 the rank-2 pertm'bation moves 
the eigenvalues of A in the horizontal direction. Since the eigenvalues i and 1 + i of A are 
on the same horizontal line, they can be simultaneously perturbed by a rank-1 perturbation. 
For t G [0,4] those eigenvalues are plotted by circles in Figure 5. We also computed the 
eigenvalues of A + tEi and of ^4 + AEi + tE2 for t € [0,4] and indicated them by asterisks 
and boxes, respectively. As is visible in Figure [5} the eigenvalues leave the straight line before 
returning to the eigenvalues of the normal matrix A + AEi + 4i?2 = A + AE. 




Figure 5: Eigenvalue trajectories of a normality preserving perturbation, and of the same 
perturbation written as the sum of non-normality preserving normal perturbations. 



6.3 Illustrations belonging to Section [5] 



Finally, we will illustrate Theorem |5. 7 The starting point is a 3 x 3 matrix A, 

and thus, UiA) = X -- 



5i 







1 

2 + 2i 




1 
2 



(84) 



With 9 = 1, the eigenvalues already lie on a parabolic curve, and thus also with the trivial 
choice Z = augmentations A+ can be constructed having eigenvalues on the same curve. 
More interesting is to choose a Z ^ such that A — iZZ* is normal. Since TiiA) has distinct 
eigenvalues, Z needs to have eigenvectors of T-LiA) as columns. Take for example 



Z 



2 

1 



then ZZ* 



5 

1 



and A — iZZ* 



i 
1 







2 + i 



The eigenvalues of A — iZZ* lie on the curve C that is the image of 



(85) 



(86) 
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Augmentations can now be constructed by choosing an arbitrary Hermitian 2x2 matrix 
M and an arbitrary unitary matrix Q, for instance 



M 



1 1 
1 2 



and Q = -V2 



and then to form the Hermitian part of as 



H{A+) 



X ZQ 
Q*Z* M 




1 


V2 
V2 





2 

■hV2 



V2 


1 
1 



(87) 



V2 


1 
2 



(88) 



Finally, A+ itself can be formed as A+ = 'H{A)+i-q{'H{A)) = 'H{A) + i-{I-'H{A))'^, resulting 



m 



A, 



5i 



V2 
V2 + V2i 




1 









2 + 2i 
hV2 



V2 



i^/2 
1 + 3^ 
l + 2^i 



V2 + V2i 


-IV2- lV2i 

1 + 2^i 

2 + 4ii 



(89) 



is normal, and computed its 



Indeed, yl+ is a 2-augmentation of We verified that A-\ 
eigenvalues. Prom this example we observe that if Z, M, Q are chosen real, then H{A^) is 
real symmetric and A^ complex symmetric, being the sum of a real symmetric matrix and 
i times a polynomial of this real symmetric matrix. Note that not all complex symmetric 
matrices are normal. In fact, the leading 4x4 principal submatrix of A^ in the above example 
is not normal, nor is the trailing 2x2 principal submatrix. Thus, A^ could not have been 
constructed using the procedure for m = 1 twice. 
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